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Cosmological Constraints on Unparticle Dark Matter 
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Abstract. In unparticle dark matter (unmatter) models the equation of state of the unmatter is given by 
p — p/(2du + 1), where du is the scaling factor. Unmatter with such equations of state would have a 
significant impact on the expansion history of the universe. Using type la supernovae (SNIa), the baryon 
acoustic oscillation (BAO) measurements and the shift parameter of the cosmic microwave background 
(CMB) to place constraints on such unmatter models we find that if only the SNIa data is used the 
constraints are weak. However, with the BAO and CMB shift parameter data added strong constraints 
can be obtained. For the TlUDM model, in which unmatter is the sole dark matter, we find that du > 60 
at 95% C.L. For comparison, in most unparticle physics models it is assumed du < 2. For the /1CUDM 
model, in which unmatter co-exists with cold dark matter, we found that the unmatter can at most make 
up a few percent of the total cosmic density if du < 10, thus it can not be the major component of dark 
matter. 



1 Introduction 

It has recently been proposed that a hidden scale- invariant 
sector of matter may exist [1] , and is named "unparticle" 
for its unusual behavior. In this scenario, there is a scale 
invariant sector with a non-trivial infrared fixed point, 
called the Banks-Zaks (BZ) field [2]. The BZ field inter- 
acts with the standard model (SM) fields via exchange of 
particles of mass Mjj : 



The unparticles do not have a fixed mass, and the den- 
sity of states of unparticle is given by 
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Typically 1 < du < 2 are considered 

ITWlTin^lTSlfn^ . For a ther- 

mal distribution of unparticles, the density and pressure 
is given by [26] 
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where Obz is the BZ operators with mass dimension dsz 
and Osm is the SM operators with mass dimension dsM- 
Dimensional transmutation produce a scale Ajj, below which 
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where Cu is the coefficient function and Ou is the unparti- 
cle operators with mass dimension du- Phenomenological 
constraints on Mu, Au, Cu, du have been derived from a 
number of particle physics [ ^ [4i r^ r6 | r7 y 8 | f W T W^ r ^ n^1 J T4 l 
[FoTfT^lfTT] and astrophysics [I^|TW2WTp2pgp4] analy- 
sis. 

Given that the unparticles interacts weakly with stan- 
dard model particles, it is natural to consider unparticle 
matter, or unmatter, as a candidate of dark matter, es- 
pecially if the unparticle could be stabilized by a discrete 
symmetry |25j . Thanks to the unusual kinematics of un- 
particles, the behavior of unmatter is distinctly different 
from the usual cold dark matter. 



where C{d v ) = B(3/2,du)r(2dv + 2)((2du + 2), 
B, _T, C are the Beta, Gamma and Zeta functions, 
equation of state for the unmatter is therefore 

wu = l/(2du + 1) 



(2) 

(3) 

and 
The 
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Thus, as the universe expands, the energy density of un- 
matter evolves as pu(z) = pua{l + Z ) 3( - 1+Wu \ If du = 1, 
this is the same as radiation, and for du —> oo its behavior 
would be similar to cold dark matter. In the intermediate 
case, its evolution would differ from both radiation and 
cold dark matter. We can then use cosmological observa- 
tions to constrain the value of du- 

We shall consider two models. In both cases we as- 
sume the universe is fiat with a cosmological constant. 
In the first case, denoted by ylUDM, the unmatter serves 
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as the sole dark matter, and the set of cosmological pa- 
rameters 9 is {J?bo, ^ua, du, ho }. In the second case, 
denoted by /1CUDM, we consider the more general case 
where the unmatter is not the only source of dark mat- 
ter. We then constrain the amount of unmatter if it does 
exist. The cosmological parameters set 9 in this case is 
{O m o, Quo, djj, ho}. We then use a Markov Monte Carlo 
Chain method to make the global fitting and constraints. 
For details of our MCMC code we refer the readers to 
Ref. [27]. 

We consider three observational constraints. The first 
one is luminosity distance moduli to type la supernovae 
(SNIa). The second is the baryon acoustic oscillation (BAO) 
feature in large scale structure as measured by the Sloan 
Digital Sky Survey (SDSS) and the Two Degree Field 
Galaxy Redshift Survey (2dFGRS). The last one is the so 
called shift parameter [28^29], which is essentially a mea- 
sure of the distance to the last scattering surface of the 
cosmic microwave background (CMB), as measured by the 
WMAP three year observation [30]. All of these provide 
constraints on the global expansion history of the universe. 



2 Methods 

The cosmic expansion rate H (z) is given by 

H 2 (z) = Hi f2(z;6), (5) 

where for AUDM: 

9) = Q w (l + zf + n AQ + f2 r0 (l + zf 

+n m {i + z? (1+Wu) (6) 

and for ACUDM: 

9) = Q mo {\ + zf + Q A q + O r0 (l + zf 

+n uo (l + zf^ +w ^; (7) 

with i?yio = 1 — ]7 m o— i7;7o — J7 r o- Here, f2 m o, f2 r o, Qao, ^uo 
are the relative abundance of matter, radiation, the cos- 
mological constant, and unmatter respectively. Of course, 
for the AUDM model, J? m0 = Qbo- 

Supernova constraint: the luminosity distance to a 
supernova is given by 

dr.M = (l + z)£*L. (8) 

and the distance moduli is 

lMh(z) = 51og 10 dz,(z) + 25, (9) 
The x 2 for the SNIa data is 

2 ,m (Hobs(zj) - Hth{Zj)f , , 

XsnW=2^ ' V 1U J 

i=i ' 

where Hobs(zi) and <ii are the observed value and the cor- 
responding error for each supernova. We use a data set of 



182 high-quality SNIa [27] selected from the Gold06 [32], 
SNLS H3 and ESSENCE [3J samples. 

CMB constraint: the CMB shift parameter R [28] 
denotes the positions of the acoustic peaks in the angular 
power spectrum of CMB, and takes the form as 

The WMAP3 data gives R = 1.70 ± 0.03 [29 , thus we 
have 

2 /R_L70\ 2 
Xr = [-QM- ) ■ (12) 

BAO constraint: we use the quantity r s /D v which is 
constrained by the BAO signature in SDSS (at z = 0.35) 
and 2dFGRS (at z = 0.2) data [3S1EE]: r s /D v (0.2) = 
0.1980 ± 0.0058 and r s /D v (0.35) = 0.1094 ± 0.0033, with 
a correlation coefficient of 0.39. Here r s is the comoving 
sound horizon size at the epoch of decoupling, and D v is 
the effective distance defined in [3T]. we do not use the 
parameter A which is extracted from the BAO measure- 
ments of the SDSS, as its definition applies to the /tCDM 
model specifically, and may not be applicable in the pres- 
ence of unmatter models [37] , 

For the combined analysis, 

X 2 =Xsn + Xr + Xbao- ( 13 ) 

We employ the Markov Chain Monte Carlo (MCMC) 
technique to calculate the posterior probability distribu- 
tions function of the parameters. The Metropolis-Hastings 
algorithm with uniform priors is used to generate the sam- 
ple, and the priors are taken as the following: .f^o 6 (0, 0.1), 
n m0 6 (0, 1), Q m 6 (0, 1), d v G (0, 10 5 ) and h £ (0.4, 0.9). 
The energy density of all components are assumed to be 
positive, fluo £ (0, 1 — fi^o/QmO — Qro) is set so that 
J?/io = 1 — Quo ^ Qbo/Qmo — Qro > 0. For each of the 
two models (ylUDM and ylCUDM) we generate six chains, 
and about ten thousands points are sampled in each chain. 
After thinning the chains, we merge them into one chain 
which consists of about 10000 points. 



3 Results 

First we consider the constraints derived purely from the 
SNIa data. For constraining the unmatter model, this is 
the most reliable one, as it is based only on the global 
expansion history, which can be calculated exactly for the 
given parameter set. 

In Fig. [T] we plot the constraint on fijj and du in the 
ilUDM model after marginalizing the other parameters. 
We found that practically all values of du are allowed. At 
large values of du, the central value of fiu is between 0.15 
and 0.25. This is what we would have expected, since for 
large value of du the behavior of the unparticle gas is very 
similar to that of the cold dark matter, and for ylCDM the 
best fit is centered in the same region. At smaller values 
of du, the contours curved to smaller fly -. 
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In Fig. [2] we plot the constraint on flu and djj in the 
ylCUDM model. Again, practically all values of du are 
allowed. For large djj where the unmatter asymptotes to 
cold dark matter, the best fit is located at flu ~ 0.22, 
as the UDM become dominant and the CDM has a very 
small abundance. However, for high values of du, all values 
of flu are allowed, in this part of the parameter space the 
UDM plays a minor role and the CDM is dominant. 
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Fig. 1. The \a, 2<r, 3<r contours of the Qu and du in the 
/lUDM model derived from SNIa observations. 
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Fig. 2. The lcr, 2<r, 3<r contours of the Qu and du in the 
ylCUDM model derived from SN la observations. 



Our constraint is improved significantly by using the 
BAO and CMB shift parameter in addition to the SNIa 
data. The r s /D v proposed by Percival et al. [3S] is a ra- 
tio of two "standard ruler", which is suitable and reli- 
able to constrain our unmatter model and could efficiently 
break the degeneracy of the parameters. The shift param- 
eter does not include all information in the CMB angu- 
lar power spectrum, but it is relatively easy to compute 



for unconventional models. The large distance to the last 
scattering surface of the CMB provides a long level arm 
for constraining the global expansion rate. We plot the 

— du contours for the ylUDM model in Fig. [31 and the 
same contours for the ylCUDM in Fig. 2J As can be seen 
from the figures, the distribution is drastically different 
from the SNIa only constraints. 

For the /LUDM model, the lowest edge of the 2a (95.5% 
C.L.) contour is about 60. Thus, small values of du, which 
is interesting from a physics perspective, is excluded at 
high confidence levels. At large du the unmatter is similar 
to cold dark matter, so again it is not surprising that the 
center value of flu is 0.23-0.29. 

For the ylCUDM model the constraint is also much 
more stringent. At 2a level, only 2 or 3 percent of the 
total cosmic density could be in unmatter. For the scal- 
ing dimension du, the minimal allowed value is 1.2 at 
95.5%C.L.. Furthermore, for the small du values, the con- 
straint on the unmatter density is stronger. For du < 2, 
the limit on unmatter density is flu < 0.01 at 2a level. 
The constraint on abundance flu loosened at greater val- 
ues of du, where the unmatter become indistinguishable 
from cold dark matter. 
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Fig. 3. The la, 2a, 3a contours of the Qu and du in the 
AUDM model derived from SNIa, BAO and CMB shift pa- 
rameter observations. 



4 Conclusion 

In this work, the scaling dimension du and the abundance 
flu of the unparticle dark matter are constrained by us- 
ing the SNIa luminosity distance moduli, BAO features in 
large scale structure as measured by SDSS and 2dFGRS, 
and the CMB shift parameter which depends on the cos- 
mic expansion history. We used the MCMC technique to 
simulate the posterior probability distributions of the pa- 
rameters. Two models are considered, viz. ylUDM (where 
unparticles are the only dark matter) and ylCUDM (where 
unparticles co-exists with other cold dark matter). 



4 



Yan Gong, Xuelei Chen: Cosmological Constraints on Unparticle Dark Matter 



tT 4 



— i — i — i — i — | — i — i — i — i — | — i — r 

ACUDM 


i i |, i i i i | 




olN laTDAUtLMB 


; / 
/ / 


- 




la 


68.3SS - 




2o 


05.5% 





3<7 




99.7% 



0.01 0.02 0.03 0.04 0.05 

Fig. 4. The la, 2a, 3a contours of the Ou and du in the 
ACUDM model derived from SNIa, BAO and CMB shift pa- 
rameter observations. 



Using only the supernova data, we find that the con- 
straints on djj and fijj are pretty weak. However, with 
the addition of the BAO and CMB shift parameter data, 
the degeneracy is broken, and strong constraint could be 
put on du and flu- For the /1UDM model, du > 60 at 2a 
level. For the ylCUDM model, if d v < 2, then Q v < 0.01 
at 95.5% C.L.; and Du is no greater than a few percent 
if du < 10. These limits severely constrained models of 
unparticle dark matter and stable relic unparticle matter. 

A major assumption adopted in this work is the equa- 
tion of state wu for the unmatter: wu = l/(2e?[/ + l). This 
expression was derived by S. L. Chen et al in Ref. [SB] for 
a thermal distribution of unparticle matter. This result is 
a little surprising, since an ideal gas of massless particles 
which also possesses scale invariance has an equation of 
state wu — 1/3. However, if one accepts the usual den- 
sity of states for the unparticles (Eq. [1]) , which was used 
for other calculations about unparticles, e.g. production 
of unparticle in particle collisions, then this equation of 
state can be derived using the standard method of statis- 
tical mechanics, as has been shown in Ref. [2B]. It is also 
possible that the unparticle matter has a non-thermal dis- 
tribution, then its equation of state would be different, and 
in that case it is not constrained by our result. However, 
if the unparticle matter was produced during the hot Big 
Bang in the standard decoupling scenario, its very nature 
must be considered a thermal distribution of unparticle 
matter. 
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